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ABSTRACT 

We examine the red fraction of central and satellite galaxies in the large zCOSMOS group catalog 
out to z ~ 0.8 correcting for both the incompleteness in stellar mass and for the less than perfect 
purities of the central and satellite samples. We show that, at all masses and at all redshifts, the 
fraction of satellite galaxies that have been quenched, i.e., are red, is systematically higher than that 
of centrals, as seen locally in the SDSS. The satellite quenching efficiency, which is the probability 
that a satellite is quenched because it is a satellite rather than a central, is, as locally, independent of 
stellar mass. Furthermore, the average value is about 0.5, which is also very similar to that seen in 
the SDSS. We also construct the mass functions of blue and red centrals and satellites and show that 
these broadly follow the predictions of the Peng et al. (2012) analysis of the SDSS groups. Together, 
these results indicate that the effect of the group environment in quenching satellite galaxies was very 
similar when the Universe was about a half its present age, as it is today. 

Key words: cosmology: observations - galaxies: groups: general - galaxies: evolution - galaxies: mass 
function - galaxies: statistics 



1. INTRODUCTION 

One of the most striking features of the galaxy popu- 
lation is the observed bimodality in the color of galax- 
ies (e.g., Strateva et al. 2001): most galaxies are either 
blue and star forming, with a star-formation rate that is 
closely linked to the existing stellar mass, or red with a 
star-formation rate that is lower by 1-2 orders of mag- 
nitude. This bimodality is observed to be already in 
place at rcdshift z ~ 1 (e.g., Bell et al. 2004; Cirasuolo 
et al. 2007; Mignoli et al. 2009) and a key question in the 
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context of galaxy evolution is to understand the mecha- 
nism^) that lead to the cessation of star-formation as a 
function of stellar mass, cosmic environment, and cosmic 
time. We refer to this process as "quenching" . 

Although the specific mechanisms of this quenching 
process are still largely debated, there has been some 
progress in phenomenologically understanding where and 
when it happens (e.g., Peng et al. 2010). There is good 
evidence that the action of quenching is different for 
galaxies in different locations with respect to the dark 
matter halos in which they reside. Broadly speaking, 
"central" galaxies are taken to be the most massive 
galaxies in a halo and to reside at the potential mini- 
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mum. "Satellite" galaxies are galaxies that are moving 
relative to the potential minimum having fallen in to the 
larger halo. In this way, all galaxies are either centrals 
or satellites. The role of the central/satellite dichotomy 
in quenching has been discussed in several observational 
papers (e.g., van den Bosch et al. 2008; Kimm et al. 2009; 
Woo et al. 2012; Peng et al. 2012). van den Bosch et al. 
(2008) introduced the "satellite quenching efficiency", 
which measures the fraction of blue (previously central) 
galaxies that are quenched when they become satellite 
galaxies. They showed that in the low-redshift Universe 
sampled by the Sloan Digital Sky Survey (SDSS), this 
quantity is independent of the stellar mass of the satel- 
lite galaxy. 

Peng et al. (2010) showed that the fraction of galax- 
ies that are quenched, i.e., the so-called "red fraction", 
in SDSS is fully separable in terms of the effects of stel- 
lar mass and environment. This means that the frac- 
tion of galaxies that survive as blue star-forming galax- 
ies can be written as the product of two functions, one of 
stellar mass and one of environment. They argued that 
this indicated the action of two distinct processes, one, 
"mass quenching", that is linked to the stellar mass of 
a galaxy but not the environment, and the other, "en- 
vironment quenching" , that is linked to its environment 
but not its mass. In Peng et al. (2012) it was shown 
that all of the environmental effects in the SDSS popu- 
lation can be ascribed to the satellite galaxies in groups. 
The colors of central galaxies arc independent of environ- 
ment, and these galaxies must experience only the mass 
quenching process. In contrast, satellite galaxies experi- 
ence both mass-quenching and environment quenching. 
It was shown in Peng et al. (2012) that the probabil- 
ity that a satellite had been environment (or satellite) 
quenched was independent of its stellar mass, as in van 
den Bosch et al. (2008), but also that it depended on the 
local density of galaxies, which was taken as a proxy for 
the location of the satellite in the group. Interestingly, 
at a given density, the satellite quenching efficiency does 
not depend on the mass of the parent halo. 

In this paper, we aim to investigate the properties 
of centrals and satellites at much higher redshifts, us- 
ing the 20k group catalog (Knobel et al. 2012, "K12") 
from the spectroscopic zCOSMOS-bright survey (Lilly 
et al. 2007). The zCOSMOS 20k group catalog contains 
about 1500 groups (including pairs) within the redshift 
range 0.1 < z < 1 and features very good and well con- 
trolled statistics in terms of completeness and purity of 
the group population. Moreover, in K12 we developed a 
probabilistic scheme to select subsamples of centrals and 
satellites which have purities up to 80%. 

The definition of the central galaxy is an important is- 
sue. In most studies, the centrals have been identified by 
just selecting the brightest or the most massive galaxies 
within dark matter halos. However, both semi-analytic 
simulations (e.g., Kitzbichler & White 2007; Henriques 
et al. 2012) as well as analysis of observational data in the 
low-redshift Universe (Skibba et al. 2011) indicate that 
a non-negligible fraction (i.e., 20%-30%) of the galaxies 
that lie at the deepest point within the potential well 
are not in fact the most luminous or most massive galax- 
ies in their halos. This ambiguity may be exacerbated 
by uncertainties in the estimation of stellar masses from 
photometric data. 



It is likely that the location of the galaxy in the group 
is more important than the property of being the most 
massive. For this reason, we developed a probabilistic 
selection scheme for centrals and satellites in K12 that 
also includes positional information in addition to stellar 
mass. This scheme was carefully calibrated using mock 
catalogs generated from numerical simulations (Kitzbich- 
ler & White 2007; Henriques et al. 2012). These well 
controlled statistics can be used, in principle, to correct 
our set of centrals an satellites for misidentifications. 

Our goal is to extend the discussion of centrals and 
satellites up to redshift z ~ 0.8 particularly with focus 
on the quenching of central and satellite galaxies as a 
function of stellar mass and cosmic epoch. We therefore 
study the fraction of centrals and satellites and construct 
the red fractions of each set. From this we derive the 
satellite quenching efficiency (averaged over different lo- 
cal densities). We also look at the mass functions of the 
blue and the red central and satellite galaxy populations 
since, as emphasized by Peng et al. (2010, 2012), this is 
also a sensitive diagnostic of the quenching action. We 
discuss our findings in the context of the model proposed 
by Peng et al. (2012). In a later study (K. Kovac et al., 
in preparation), we will analyze the quenching of cen- 
trals and satellites in zCOSMOS with emphasis on their 
dependence on the local environmental density. 

The layout of this paper is as follow. In Section 2, we 
describe the data and the samples that we used for our 
analysis. In Section, 3 we describe our results for the 
fraction of satellites, the satellite quenching efficiency 
and the mass functions of the different galaxy popula- 
tions. Finally, in Section 4 we summarize and conclude 
our findings. 

Where necessary, a concordance cosmology with Hq = 
70 km s -1 Mpc" 1 , fi m = 0.25, and fi A = 0.75 is applied. 
All magnitudes are quoted in the AB system. We use 
the term "dex" to express the antilogarithm, i.e., 0.1 dex 
corresponds to a factor 10 ' 1 ~ 1.259. By log we refer to 
the 10-based logarithm and we occasionally use log M as 
short term for log(M/AfQ). 

2. DATA 

We use the central and satellite classification scheme 
from the group catalog described in K12, which is based 
on the 15 < Iab < 22.5 flux limited sample from the 
spectroscopic zCOSMOS-bright survey (Lilly et al. 2007, 
2009, Lilly et al. 2012 in preparation). Multiplying the 
(random) spatial selection with the redshift success rate, 
this sample has an overall mean sampling rate of the tar- 
get sample of about 50%. The catalog was produced by a 
Friends-of-Friends (FOF) multi-run scheme (see Knobel 
et al. 2009) in which successively different group-finding 
parameters were used, optimized for different richness 
groups. For details on the construction and the statis- 
tics of the group catalog we refer to K12. 

As well as attempting to classify all of the galaxies 
as centrals or satellites, taking into account also those 
galaxies for which a spectroscopic redshift is not avail- 
able, a feature of this group catalog is the possibility 
to also produce fairly pure subsamples of centrals and 
satellites. The purity of a (sub)sample is defined as the 
fraction of correctly classified galaxies (either centrals or 
satellites) within the (sub)samplc divided by the total 
number of galaxies in the (sub)samplc. In K12 every 
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galaxy (including those only with photometric redshifts, 
photo-z) was assigned a membership probability p to be 
a member of a given group, a probability pm to be the 
most massive galaxy within that group, and another pa- 
rameter pma which additionally to pm also considers the 
positional information within the group. The photo-z 
galaxies were used in K12 to help classify the spectro- 
scopic galaxies, e.g., by indicating the possible presence 
of a better candidate for the central, but are not then 
used in the present analysis because their own purities 
are naturally much lower. 

By selecting (spectroscopic) galaxies according to their 
p, pm, and pma, (sub) samples of centrals and satellites of 
varying size and purity can be produced. An important 
case is the "dichotomous sample" , in which all galaxies of 
the zCOSMOS spectroscopic sample are classified to be 
either a central or a satellite galaxy. Based on simulated 
mock catalogs, it was found that defining the satellites by 
simultaneously applying p > 0.1, pu < 0.5, and pua < 
0.5, and identifying all other galaxies (including those 
that are not associated to any group, i.e., with p = 0) 
as centrals, resulted in a purity (averaged over the whole 
population) of about 84% for the centrals and about 74% 
for the satellites (see Table 10 of K12). 

It should be appreciated that these purities include 
both the difficulties of identifying the central of a given 
group, and the errors that are made in identifying the 
groups in the first place. Clearly, a group that has been 
fragmented into multiple structures, or over-merged with 
other groups into a spuriously large structure by the 
group-finder algorithm (see K12 for discussion), will not 
end up with a 100% correct classification of centrals and 
satellites. 

The impurity of the adopted central and satellite pop- 
ulations due to cross-mixing between them, will lead to a 
reduction of any observed differences in a given property 
between the two sets. If the purities are known, then of 
course this dilution effect can be corrected for, as detailed 
below. We would expect the purities to be a function of 
both redshift and stellar mass, and therefore computed 
them in this way. The purity functions for the dichoto- 
mous sample are shown in Figure 1. For both centrals 
and satellites the purity is > 80% at all redshifts and is 
essentially independent of mass. 

It should be noted that the purities cited in K12 were 
obtained using mock catalogs that were derived from the 
light cones of Kitzbichler & White (2007), which are 
in turn based on the Millennium Simulation (Springel 
et al. 2005) and the semi-analytic model of Croton et al. 
(2006). For our present analysis we updated our mock 
catalogs using the light cones of Henriques et al. (2012), 
which are also derived from the Millennium Simulation, 
but are based on the more recent semi-analytic model of 
Guo et al. (2011). The differences on the derived purities, 
however, are rather small and our final results are essen- 
tially independent of the specific semi-analytic model, 
from which our mock catalogs were derived. 

2.1. Data samples 

We perform our analysis in three independent redshift 
bins: 0.1 < z < 0.4, 0.4 < z < 0.6, and 0.6 < z < 0.8 
(see Table 1). In each redshift bin, the galaxy sample 
is divided into blue (assumed star forming galaxies) and 
red (assumed quenched) galaxies. We apply the same 



Table 1 

Galaxy samples 



Redshift 


M ■ a 
iw min 


# centrals 


# satellites 


red blue 


red blue 


0.1 < z < 0.4 


9.9 


788 788 


343 228 


0.4 < 2 < 0.6 


10.1 


746 574 


187 83 


0.6 < z < 0.8 


10.5 


1122 359 


212 53 



Note. — The numbers refer to the dichotomous sam- 
ple (see the text) and consider only galaxies with masses 
larger than M m ; n . 

a Minimum stellar mass in units of log(M m j n / Mg) . 



color cut to all three redshift bins, where red galaxies 
are defined by 

U-B> l.l + O.O2(log(M/M )-lO) (1) 

with U and B being absolute magnitudes and M being 
the stellar mass. This cut was obtained by eye from the 
overall color-mass diagram at each redshift. Changing 
the cuts by 0.02 in color, or making small changes to the 
slope, have small effects on the results, typically much 
less than the quoted la error bars. 

The absolute magnitudes and stellar masses were de- 
rived from spectral energy distribution (SED) fitting us- 
ing ZEBRA+ (sec K12 for details). The masses arc defined 
here to be the integral of the star-formation history of the 
model galaxy, and do not, therefore, include the effects of 
mass return from the stellar population. The advantage 
of this choice is, for instance, that red (passive) galaxies 
with zero star formation rate have constant mass, en- 
abling straightforward comparisons of stellar mass func- 
tions. 

Since our sample is /-band flux limited, the complete- 
ness of our sample depends on stellar mass and on the 
redshift. Moreover, the mass-completeness is different 
for red and blue galaxies. Given the excellent photo-z 
in the COSMOS field, we use the photo- z sample down 
to 7ab < 24 to empirically determine and correct for 
this incompleteness. We apply the same color cut to all 
COSMOS galaxies and compute for each redshift bin, 
and for both red and blue galaxies separately, the com- 
pleteness as a function of stellar mass, i.e., what fraction 
of galaxies in the photo- z sample are also present in the 
spectroscopic sample. To compute fractions of galaxies 
and mass functions, individual objects in the spectro- 
scopic sample are then weighted by the inverse of the 
completeness. The mass-limit M m i„ of the survey at a 
given redshift (see Table 1) is defined to be that mass at 
which the weighting reaches a factor of two relative to 
the high mass end. We checked that our results are not 
sensitive to the exact choice of M mm . 

This correction method is based on the assumption 
that, within a given color class (blue or red) the mass-to- 
light ratio of a galaxy does not depend on whether it is 
a central or satellite, i.e., on whether it is a group mem- 
ber or not, so that the correction can be applied in the 
same way to centrals and satellites. Clearly the color- 
mix between red and blue is different for centrals and 
satellites, but we know of no evidence for such variation 
within each color-class. For instance, Peng et al. (2012) 
showed that the specific star- formation rates of blue star- 
forming satellites and of blue star-forming centrals were 
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essentially indistinguishable, even though their numbers, 
relative to passive galaxies, varies markedly. 

Statistical errors are computed by bootstrapping. 
That is, to compute the error for a certain quantity (e.g. a 
fraction of galaxies or a mass function), we randomly re- 
sample our whole galaxy sample 100 times, perform the 
whole analysis consistently for each of these samples, and 
finally compute the dispersion (e.g., standard deviation 
or quartiles) of the corresponding results. 

3. RESULTS 

In this section, we develop the results from our anal- 
ysis. To make the formalism as simple as possible, we 
adopt the following conventions. For any quantity the 
subscript 'r' refers to 'red', 'b' to blue, 'c' to central, 
and 's' to satellite. If subscripts are combined, they re- 
fer to the intersection of the corresponding subsamples, 
i.e., / r x denotes the fraction of red centrals. Whenever 
a quantity is given a tilde, it is uncorrected with respect 
to color completeness or central/satellite impurities. 

3.1. Fraction of satellites 

Although the purity for centrals and satellites is quite 
good for all samples (see Figure 1), the impurities are 
still on the level of 20% and this will affect our results. 
For this reason, we will use these purity functions to 
correct all of our results, as well as possible, for this 
contamination. As an example, the dichotomous galaxy 
sample will consist of N c galaxies classified as centrals, 
and N s as satellites, within a given bin in stellar mass 
and redshift. If P c is the purity of centrals and P s is the 
purity of satellites within that bin, the corrected numbers 
of centrals and satellites are given by 

N c = N c - (1 - P C )N C + (1 - P S )N S (2) 

N s = N s - (1 - P S )N S + (1 - P C )N C , (3) 

respectively. 

The fraction of satellites as a function of mass is shown 
in Figure 1 in the three redshift bins and is compared 
to the mock catalogs that were derived from the sim- 
ulated light cones of Henriques et al. (2012). The red 
lines show the purity-corrected satellite fraction for our 
actual zCOSMOS data, the black lines show the purity 
corrected satellite fraction for the corresponding sample 
in the mock catalogs, and the gray shaded area shows the 
la range for the real fraction of satellites in the mock cat- 
alogs. The dashed lines show the fraction of satellites for 
the mock sample, if we do not correct for impurities, and 
show the size of effect that is being corrected for. It is 
clear that the satellite fraction would be systematically 
too low compared to the correct result (gray area) if the 
purity correction is not applied. The fraction of satel- 
lites for the actual data, once corrected, is in excellent 
agreement with the expectations from the cosmological 
simulations. We find a satellite fraction between 20%- 
40% in all redshift bins. 

Comparing the result in the three redshift bins, we de- 
tect a small, but significant decrease in the satellite frac- 
tion f s with redshift at the lower masses (a similar result 
integrated over all masses was shown in K12). For the 
lower masses (i.e., M < 1O 1O 5 M0), / s increases system- 
atically from z ~ 0.7 to z ~ 0.3 by up to a factor of two, 
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Figure 1. Fraction of satellites / s as a function of stellar mass M 
within the three redshift bins. The red lines show the (corrected) 
fraction for the zCOSMOS data with the error bars being the stan- 
dard deviation derived from bootstrapping. The black lines show 
the fractions for the corresponding mock samples selected in the 
same way as the actual data, where the solid lines correspond to 
the fractions corrected for impurities (see the text) and the dashed 
lines to the uncorrected fractions. The error bars represent the 
standard deviation of the 24 mock catalogs. The gray shaded areas 
represent the lcr regions of the fraction of real satellites within the 
24 mock catalogs. The blue solid lines show the purity of centrals, 
and the blue dashed lines the purity of satellites. The agreement 
between the actual data and the mock catalogs is very good. 

while at the high mass end there is no significant evo- 
lution detected. The increase in the satellite fraction is 
as expected from the hierarchical buildup of groups over 
cosmic time. Comparing our results with that of van den 
Bosch et al. (2008) from SDSS, we find that the satellite 
fraction observed by them is systematically lower, by at 
least 25%, than the satellite fraction in our lowest red- 
shift bin and decreases more strongly with mass. These 
differences could be at least partially explained by the 
fact that van den Bosch et al. (2008) did not correct for 
impurities in their sample and thus might underestimate 
the real fraction of satellites (cf. the dashed line in Figure 
1)- 

3.2. Satellite quenching efficiency 

It has been shown in the low-redshift Universe that 
the red fraction of the centrals / rjC (M) at a given stel- 
lar mass M is lower than the red fraction of the corre- 
sponding satellite galaxy population / rjS (M) (e.g., van 
den Bosch et al. 2008; Peng et al. 2012). It should be 
appreciated that in comparing centrals and satellites at 
the same stellar mass, the "centrals" in question will gen- 
erally be in different structures (i.e., lower mass halos) 
than those which contain the "satellites" , i.e., the central 
galaxies of the halos containing these particular satellites 
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will generally be of higher stellar mass. However, our in- 
terest is in comparing galaxies of the same stellar mass, 
so as to remove the (strong) effects of mass- quenching. 

Using the red fractions of central and satellite galaxies, 
we can straightforwardly compute the average "satellite 
quenching efficiency" e s defined as 



3 (M) 



/ r , B (M)-/ r , c (M) 
A,c(M) 



(4) 



The average in this case is taken over all environments 
within the halos, i.e., the local densities or group-centric 
distances. e s can be interpreted as the fraction of (sur- 
viving) blue centrals at a given stellar mass that are 
quenched when they become satellites by falling into an- 
other dark matter halo. The value as measured in SDSS 
is about 0.4 (uncorrected for purity) and is essentially 
constant with mass (van den Bosch et al. 2008; Peng 
et al. 2012). However, as mentioned above, the precise 
value may depend on the specific definition of central 
and satellite galaxies, and whether one corrects for im- 
purities. 

In order to compute / rjC (M) and f T>s (M) in the zCOS- 

MOS sample, we again correct for impurities. If / r . c and 

/r, s are the observed (uncorrected) red fractions of cen- 
trals and satellites, in a given stellar mass bin, then these 
will be related to the true red fractions / r c and / rjS of 
centrals and satellites by 



/r,c /r.c^c ~t~ (1 -fc)/r,s: 
/r,B - /,.,/'. + (1 - Ps)fr,c- 



(5) 
(6) 



These simultaneous equations can be trivially solved to 
give 



/r,s = ^^/r,c(l ^s) fr,sPcJ 



where 



C = (1 - P c )(l - P s ) - P C P S . 



(7) 

(8) 

(9) 

Clearly / r c and / rjS are not well defined if both P c and 
P s are equal to 50%, as C is zero in this case, and there 
is no way to recover the correct red fractions if both the 
sample of centrals and the sample of satellites are equal 
mixtures of (real) central and (real) satellites. It should 
be noted that the above equations are general and may 
be applied for subsets of centrals and satellites, chosen to 
maximize their purities, and not just to the dichotomous 
sample in which every galaxy is retained, provided that 
the construction of the subsample(s) is not dependent on 
the color of the galaxy. From this point on, all quoted 
red fractions are corrected in this way. 

The (corrected) red fractions of centrals and satellites 
for the dichotomous sample are shown in Figure 2 for all 
three redshift bins. The red fractions of the satellites / rjS 
(blue lines) are always higher than those of the centrals, 
at all redshifts and at all masses probed by this study. 

The corresponding satellite quenching efficiency e s 
(black solid lines) is also found to be essentially indepen- 
dent of mass, at each redshift, and to scatter around a 
value of ~ 0.5. Within our sample there is no clear change 
in this quantity with redshift. Furthermore, the value is 
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Figure 2. Red fraction of centrals and satellites as a function of 
stellar mass M within the three redshift bins. The red lines re- 
fer to the centrals and the blue lines to the satellites, respectively, 
of a dichotomous sample. The solid black line corresponds to the 
satellite quenching efficiency derived from these red fractions. For 
comparison, the dashed black lines show the satellite quenching 
efficiencies for high purity samples of centrals and satellites. At 
any redshift the red fractions of centrals and satellites are clearly 
distinct and the satellite quenching efficiencies are essentially in- 
dependent of mass. All error bars represent the upper and lower 
quartiles derived from bootstrapping within each mass bin. All 
quantities are corrected for color incompleteness and impurities 
due to misidentifications of centrals and satellites. 



consistent with that found at low redshift in SDSS (van 
den Bosch et al. 2008: Peng et al. 2012) especially if the 
SDSS measurements are corrected (as they should be) for 
impurities in the samples. We note that there may be ev- 
idence for a possible deviation from mass-independence 
at the very highest stellar masses, where the red fraction 
of satellites in two of our three redshift bins increases 
sharply towards unity, but the statistical significance of 
this effect is not yet compelling. It may also be present 
in the SDSS analyses of van den Bosch et al. (2008) and 
Peng et al. (2012), although the error bars are also big 
in these analyses. 

To test our results for robustness, we recomputed e s 
for the (corrected) red fractions of high purity samples 
of centrals and satellites. A high purity subsample of 
centrals is just given by the set of spectroscopic galaxies 
that are not associated to any group (i.e., with member- 
ship probability p = 0), and a high purity subsample of 
satellites is given by selecting all galaxies with p > 0.5, 
Pm < 0.1, pma < 0.1 that are in groups with at least four 
observed spectroscopic members. Both of these samples 
are about 2%-7% purer than the corresponding dichoto- 
mous samples, where the difference in purity increases 
with mass. The gain in purity comes at a cost in com- 
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pleteness, and these samples are about 20%-50% smaller 
than the dichotomous sample depending on mass and 
rcdshift. The resulting e s (dashed black lines) are quite 
similar to those from the dichotomous sample, within the 
statistical errors, and thus support the robustness of our 
results. Only in the lowest mass bin at high rcdshift is 
the difference of order of 2a. We also checked that the 
mass- independence of e s is not sensitive to our correc- 
tions (either for the color completeness or for impurities), 
as the corresponding e s derived from the uncorrected red 
fractions is just systematically shifted down, albeit by 
about 60%. 

Thus we conclude that our results are consistent with 
an essentially mass and redshift independent satellite 
quenching efficiency e s , with a value of about 0.5 (with a 
scatter of this value from redshift bin to rcdshift bin of 
< 0.1), and that e s is not substantially affected by sys- 
tcmatics. It should also be noted that our mean value 
of 0.5 may well be consistent with the measurements in 
the low-redshift Universe that have yielded a measured 
value of 0.4, if impurities of the order 10% are present 
in the corresponding low-redshift samples of centrals and 
satellites. 

The current analysis therefore extends the previous re- 
sults about the (average) quenching of satellites, reported 
in SDSS, back to the epoch when the Universe was about 
a half of its present age. It is quite striking that the prob- 
ability of a satellite to have been quenched (because it 
is a satellite) is the same over a wide range of cosmic 
epochs, in addition to it apparently being independent 
of the stellar mass. This must provide a clue as to the 
physical nature of satellite-quenching. 

3.3. Stellar mass function 

As pointed out in Peng et al. (2010) the stellar mass 
functions of the red and blue galaxy population can give 
a clear view of quenching processes. Not least, Peng et 
al's mass-quenching process is the process which deter- 
mines the characteristic mass M* of the (surviving) star- 
forming population, and establishes the M* and a of the 
passive mass-quenched galaxies in terms of the values of 
the star-forming population. 

In this section, we therefore compute the mass func- 
tions for centrals and satellites split into red and blue 
galaxies. In deriving the mass functions, we use the same 
volume for centrals and satellites, and so the mass func- 
tions of the two samples will straightforwardly add to the 
overall mass function of galaxies in the Universe. 

To compute the mass functions for the different sub- 
populations we again correct for the mass-incompleteness 
in the two color bins, and for the central/satellite- 
impurities. Similar to Equations (2) and (3), we correct 
for the impurities as follows 
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where the red (and accordingly the blue) fractions are 
given by Equations (7) and (8). The corrected mass func- 
tions for the centrals and satellite galaxies of the red (red 
points) and blue populations (blue points) in the three 
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Note. — The values correspond to the best fitting param- 
eters of a single Schechter function (14) to the mass functions 
in Figure 3, and the error bars indicate their estimated 68% 
confidence regions. We want to point out that the error bars 
of the three parameters for a given fit are strongly correlated 
(cf. Figure 4). 

rcdshift bins are shown in Figure 3. All mass functions 
can be well fitted by Schechter-functions 



dlogM 



M 
~M* 



cxp 



(14) 



where In = log e denotes the logarithm to the basis e, 
and the best fitting parameters (corresponding to the 
dashed curves in Figure 3 and the crosses in Figure 4) 
are summarized in Table 2. 

The ler confidence intervals for the Schechter fits in the 
M*-a-plane are shown in Figure 4, where we projected 
out the (j)* parameter, since for this analysis we are not 
interested in the normalization of the mass functions. We 
clearly see the following trends: First, the error contours 
increase in area with redshift as the available range of 
mass shrinks. Second, the parameters of the red compo- 
nents of both centrals and satellites are typically better 
constrained than the parameters of the blue mass func- 
tions due to the larger numbers of centrals in all samples 
and the greater curvature of the mass functions due to 
the less negative a. Third, the mass functions of blue 
satellites are particularly poorly constrained. In fact, at 
medium and high redshift the information on the blue 
satellite mass functions is quite limited. 

It is well established that the stellar mass function of 
the overall galaxy population in the low-redshift Universe 
is a double Schechter function (Baldry et al. 2006, 2008; 
Peng et al. 2010, 2012; Baldry et al. 2012) and this also 
extends out to z ~ 1 (Pozzetti et al. 2010; Bolzonella 
et al. 2010). Peng et al. (2010, 2012) provided a natural 
explanation for this in terms of the action of their mass- 
quenching. The fact that the M* (and a) of star-forming 
galaxies is observed to be constant to at least z ~ 2 (e.g., 
Ilbert et al. 2010; Peng et al. 2010), despite the very large 
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Figure 3. Stellar mass functions for the red and blue populations of centrals and satellites within the three redshift bins. The red points 
correspond to the red populations, the blue points to the blue. The dashed lines show the best Schechter fits (see Table 2) and the solid 
lines the best Schechter fits for a fixed M* = ifj 10 - 95 Mq (see the text). The black points and lines correspond to the sum of the blue and 
red components. All error bars represent the standard deviation derived from bootstrapping. 



0.1<z<0.4 0.4<z<0.6 0.6<z<0.8 




-2 -10-10-10 

a a a 

Figure 4. lcr confidence regions in the a-M* plane from the 
Schechter fits to the mass functions. The third Schechter param- 
eter (ft* is projected out. The red contours correspond to the red 
galaxy population and the blue contours to the blue one. The solid 
contours refer to centrals and the dashed ones to satellites. The 
crosses mark the best fitting parameters in the plane. All con- 
fidence regions at all redshifts are essentially consistent with the 
same M* ~ m 10 - 95 Mq indicated by the dashed line. 



increase in stellar mass implied for a galaxy that remains 
on the star-forming Main Sequence, places a very strong 
constraint on the form of mass-quenching (sec Peng et al. 
2010 for details). This in turn produces well-defined mass 
functions of the passive galaxies. 

The approach of Peng et al. (2010) predicts precise 
quantitative relationships between the Schechter param- 
eters of different components of the galaxy population. 
Applied to the central/satellite population, as in Peng 
et al. (2012), these become: 

• The mass functions of blue centrals and blue satel- 
lites should be single Schechter functions with the 
same M* and a and differing only by their normal- 
ization <f>* . 

• The mass function of red (passive) centrals is pro- 
duced by mass quenching alone: It should be a 
single Schechter function with the same M* as the 
star-forming galaxies, but with a faint end slope 
that differs in a by about unity. 

• Both mass quenching and environment quenching 
act on satellite galaxies. As a result, the mass func- 
tion of red (passive) satellites should be a double 
Schechter function with the same M* , but with 
two different slopes, the one with a + 1 that is pro- 
duced by mass quenching (as for the centrals) and 
another with the same a as the star-forming popu- 
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lation that is produced by the environment quench- 
ing, which recall is independent of mass. 

Thus, all four mass functions are characterized by the 
same M* and by a faint-end slope that is either a or 
a + 1 (or a superposition of them in the case of the red 
satellites). Merging of galaxies can cause small modifi- 
cations of these relations, but these are evidently small 
and confined to a small increase in M* of passive central 
galaxies and, in the main, these relations are observed 
with impressive precision (Peng et al. 2012). 

It is of clear interest to see whether these relations 
still hold for our set of centrals and satellites at redshifts 
extending out to unity. We first simply fit all of the 
components with a single Schechter function. It can be 
seen in Figure 4 that all the Schechter components are 
consistent, at all redshifts, with a single value of M* ~ 
1O 1O ' 95 M0 (dashed line), which in turn is consistent with 
the values of M* for star forming galaxies observed in 
COSMOS (see Peng et al. 2010, Table 1). It should be 
noted that the Peng et al. (2010, 2012) values for the 
SDSS mass functions quote a ^0.2 dex lower value of M* 
because they were based on calculations of stellar masses 
that included the effects of mass return from stars. 

Since M* is degenerate with a, we test the predictions 
for a from the model of Peng et al. (2012) by fixing M* 
at the value of 10 10 95 Mq, since the model predicts the 
same M* for all mass functions. The resulting values of 
a are shown in Figure 5 and the corresponding Schechter 
models are shown in Figure 3 as solid lines. We see a clear 
difference between the faint-end slope of the red popula- 
tion and the faint-end slope of the blue population, and 
the difference Act is close to unity in each case at all red- 
shifts. Unfortunately, in Figure 3 we cannot well isolate 
the second Schechter component for the red satellites, as 
it becomes most apparent at M < 10 10 Mq in SDSS. 
However, the fact that satellites have a different color 
distribution to centrals (see Figure 2) implies that the 
effect of the secondary, satellite-quenched, population of 
passive satellites should be present. We therefore will 
assume a double Schechter function for the passive satel- 
lites whose components have the same M* = 1O 1O 95 M0, 
have a differing by unity, and have a relative </>* that is 
set by Equation (11) of Peng et al. (2012). This is shown 
as the second points (open triangles) for the passive satel- 
lites in Figure 5. 

Fitting then for just the a of each passive population, 
we find a ~ —0.3 for red centrals and a slightly lower 
a for red satellites at any rcdshift. In each case, the 
blue population mirrors the passive one with a faint end 
slope a that is steeper by about unity (i.e., Aa ~ 1). In 
the model of Peng et al. (2012) the red mass functions 
are set from the blue ones. Here, we consider the other 
way around just because the red faint end slopes are 
better determined observationally. The fact that the a 
for the satellites is always a bit smaller than the a for 
the centrals reflects the gradual decrease of the satellite 
fraction / a with stellar mass of the galaxy (cf. Figure 1, 
see Section 6.2 of Peng et al. 2012 for a discussion of 
this effect). We checked that our values for a are not 
sensitive to the specific values of M m j n (see Table 1), i.e., 
changing M m - ln by ±0.2 dex leads to a shift in Aa that 
is typically of the order of the error bar. 

We conclude that in all of our redshift bins, the ob- 
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Figure 5. Faint-end slope a if M* is fixed at 10 10 ' 95 Mq. The 
squares correspond to centrals and the triangles to satellites, where 
each time the red symbols correspond to the red galaxy population 
and the blue symbols to the blue one. The filled symbols refer to 
single Schechter function fits, while the open triangles correspond 
to double Schechter function fits (see the text). The error bars 
indicate the standard deviation. For both centrals and satellites 
the differences in the faint-end slope Aa between the red and the 
blue components are about unity at any redshift. 

served trends are very similar to those described in Peng 
et al. (2012) as observed with SDSS (see their Figure 13 
and Table 2), within the obvious limitations of the data, 
and offer further support for the validity of their sim- 
ple phenomenological model of the galaxy population, 
and specifically the role of centrals and satellites out to 
z ~ 0.8. 

4. CONCLUSIONS 

The goal of this paper was to examine the population 
of central and satellite galaxies at redshifts out to z ~ 0.8 
as seen in the extensive group catalog that we have earlier 
constructed from the zCOSMOS-bright rcdshift survey. 
Specifically we have looked at the fraction of centrals and 
satellites that have been quenched (i.e., in which star- 
formation has effectively ceased) in order to see if the 
central/satellite paradigm that has been established in 
the local Universe by studies of the SDSS is maintained 
at these significantly earlier epochs. Our analysis has 
been framed by the SDSS analysis of Peng et al. (2012), 
which itself was developed from the quenching formalism 
of Peng et al. (2010), and the reader is referred to those 
papers for detailed discussion. 

Our analysis has carefully corrected for both the differ- 
ent mass incompletenesses of blue and red galaxies (cal- 
ibrated from photo- z data) and for the impurities in the 
samples of centrals and satellites that enter due to both 
ambiguities in the identification of the central galaxies 
and through imperfections in the input group catalog. 
These impurities are calibrated from comparison with 
mock samples. For the analysis of the red fractions we 
are able to work with both a dichotomous sample, in 
which all spectroscopic galaxies are classified as cither 
central or satellite, and with subsamples of centrals and 
satellites that are expected to have higher purity. 
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We find that the satellite quenching efficiency, which 
is the probability that a galaxy is quenched because it 
is a satellite, and which is obtained by comparing the 
red fractions of centrals and satellites at the same stellar 
mass, is, as locally, independent of stellar mass (except 
possibly at the very highest masses). Furthermore, the 
strength of the satellite quenching efficiency, averaged 
over all satellites, is very similar in zCOSMOS, out to 
z ~ 0.8, as in SDSS locally, at a value of about 0.5. 

We also find that the mass functions of red and blue 
central and satellite galaxies broadly follow the expecta- 
tions established by the Peng et al. (2010, 2012) model, 
i.e., within the uncertainties they are consistent with the 
same characteristic M* and the correct relations between 
the faint-end slopes a of the various populations. This 
confirms the applicability of a universal mass-quenching 
process applying to both satellite and central galaxies. 

Taken together, these results go a long way to confirm- 
ing the operation of the different quenching channels at 
z ~ 0.8, i.e., the mass-quenching process applying to all 
galaxies, and the satellite- (or environment-) quenching 
applying to just satellite galaxies. Not least it is quite 
striking that the probability for a satellite galaxy to be 
quenched (relative to the case if it were a central) given 
by the (density-averaged) satellite quenching efficiency e s 
is apparently unchanged at z ~ 0.8 as today. Of course 
the corrections applied to derive these results are signif- 
icant and precise quantitative values should be treated 
with caution. 
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